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We demonstrate the possibility to efficiently split the near-field heat flux exchanged between 
graphene nano-disks by tuning their doping. This result paves the way for the developement of an 
active control of propagation directions for heat fluxes exchanged in near-field throughout integrated 
nanostructures networks. 


The control of electric currents in solids is at the origin 
of modern computer technology which has revolutionized 
our daily life. Until the 2000s no thermal counterpart 
had been developed to control the flow of heat at the 
nanoscale in a similar manner. In 2006, a step forward 
in this direction has been done by Li et al. [1] when they 
introduced the first concept of a thermal transistor for 
controlling heat fluxes carried by phonons through solid 
segments. Later, several prototypes of phononic ther¬ 
mal logic gates [2] as well as thermal memories [3] were 
developed in order to process information by means of 
the heat fluxes carried by phonons [3j. Besides, different 
solid-state thermal diodes were conceived nm allowing 
for rectifying these fluxes in asymmetric solid segments. 

Very recently, there has been a fast growing interest 
in developing active functionalties to manage heat trans¬ 
fers by radiation rather than by conduction between con¬ 
tactless solids. Since 2010 several radiative thermal rec¬ 
tifiers hohh] have been proposed using spectrally selec¬ 
tive nanostructures and phase-changes materials. In 2014 
and 2015 a radiative analog of a transistor was suggested 
theoretically which allows for switching, modulating and 
amplifying the heat flux exchanged both in the near¬ 
field [2T] or far-field regime (22]. Furthermore, a concept 
of a radiative thermal memory was introduced working 
in the far-field |23J and the near-field regime [24] . Such 
devices open the way for new perspectives concerning the 
development of contactless thermal circuits intended for 
an active thermal management with photons instead of 
electrons or phonons. Finally, the thermal diode con¬ 
cept based on phase-change materials has already been 
tested, experimentally, in the far-field regime [20] in 2014. 
A review of these recent developments can be found in 
Ref. [25] . 

In this Letter, we introduce the concept of a heat flux 
splitter which allows us to tune the direction of propaga¬ 
tion of heat flux exchanged in the near-field. To demon¬ 
strate the operating modes of this heat splitter we con¬ 


sider a set of three graphene nano-disks in mutual inter¬ 
action. We show that the heat flux exchanged between 
these nano-disks cannot only be splitted equally in two 
predefined directions as for a 50:50 beam splitter, which 
is trivial, but it can also be oriented in mostly one pre¬ 
defined direction acting like a 99:1 beam splitter by an 
appropriate tuning of the Fermi levels in the graphene 
nano-disks. Since this tuning can be achieved by electri¬ 
cal gating, for instance, we thus demonstrate, in particu¬ 
lar, that the direction of propagation of the radiative heat 
flow can be dynamically controlled in nano-architectures 
by electrical means. 
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Figure 1: Sketch of the heat flux splitter. Three graphene 
disks with different Fermi levels exchange thermal energy in 
the near-held through many-body interactions. The magni¬ 
tude of heat flow can be controlled by an appropriate tuning 
of graphene Fermi level. 

The geometrical setup of the heat flux splitter is 
sketched in Fig. 1. Three nano-disks of graphene are 
distributed inside a transparent dielectric of permittivity 
eh • Each nano-disk is located at its mass center posi¬ 
tion Vi and its orientation with respect to the canonical 





2 


frame ( x , z^) which is determined by the rotation 
matrix R^. Those nano-disks are supposed to be held 
at different temperatures Zb Furthermore, for the sake 
of clarity, we assume here that the nano-disks are small 
enough compared with the smallest thermal wavelength 
A Ti = ch/(k&Ti) (c is the vacuum light velocity, 27 rh is 
Planck’s constant, and k b is Boltzmann’s constant) so 
that they can be modeled by simple radiating electrical 
dipoles. The polarizability of graphene nano-disks in the 
and directions is described by a Lorentzian line 
shape [26] 
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where uo p ^ is the plasma frequency of the i th disk, Ki is 
its decay rate, and K r ^ is the radiative contribution to 
Ki. Along the z^ direction this polarizabiliy vanishes 
so that, in the canonical basis, the polarizability tensor 
reads = diag(c^, a$, 0). Thus, the polarizability ten¬ 
sor for the oriented nano-disk is simply described in the 
canonical frame by the matricial product 

a,( W ) = R,a/°)( W )R- 1 , (2) 

where denotes the rotation matrix related to the i th 
disk (see Fig. 1). 

The radiative power exchanged between the nano-disks 
can be described using the many-body radiative heat 
transfer theory 1271 Eg extended to anisotropic struc¬ 
tures El Eg. By neglecting the far-held radiative losses 
towards the surrounding, which is typically small com¬ 
pared to the near-held contribution, the net power ex¬ 
changed in the near-held between the i th and the j th 
nano-disk reads m 

r oo i 

Vij = 3 J £ [0(W, Ti) - Q(co, Tj)]Tij(co), (3) 

huj 

where @(u;,T) = fuj/[e k B T — 1 ] is the mean energy of a 
harmonic oscillator in thermal equilibrium. denotes 
the transmission coefficient between the i th and the j th 
nano-disk and is dehned as 


%,j (w) = 2ImTr , (4) 

where Xj: &ij and are the susceptibility tensor plus 
two matrices which are given in terms of Green’s ten¬ 
sors and the polarizability tensors. These quantities are 
dehned in Ref. [29] . 

Let us now apply this theoretical framework to describe 
the operating modes of near-held heat hux splitter. To 
this end, we consider a triplet of equal graphene disks 
of diameter D as depicted in the inset of Fig. 2. The 
center of mass of the hrst disk is located in the origin 
of the canonical frame while the centers of mass of the 
two other disks are located on a virtual sphere of radius 




Figure 2: (a) Thermal power exchanged in the near-held be¬ 
tween graphene disks lOOnm radius versus the separation dis¬ 
tance in a three body system. The Fermi level of disks 1, 2 
and 3 are Ef x — 0.1 eV, Ef 2 = 0.1 eV and Ef 3 = 0.8 eV while 
their temperatures are Ti = 350 K and T 2 = T 3 = 273 K. All 
disks are parallel to the x^-y^ plane and 6h — 1. (b) Plot 
of V 12 KV 12 + V 13 ) and V 13 KV 12 + V 13 ) for the same config¬ 
uration. 


R which is resting on the plane such that its 

contact point coincides with the origin. For convenience 
we consider only the situation, where all centers of masses 
are bound to the x^-z^ plane and equi-distant to each 
other. Hence, the distance between the mass centers of 
the disks is l = ^-R. 

Now, let us hrst consider the case, where all disks are 
parallel to the x^-y^ plane. Furthermore, we assume 
that T\ = 350 K and Z 2 = Z 3 = 273 K so that there is 
a heat hux from disk 1 (see Fig. 2) to disk 2 and 3. If 
the Fermi levels of all three disks are the same then the 
power V 12 an d Z 13 exchanged between disk 1 and 2 and 
between 1 and 3, respectively, will be the same as well. 
This situation corresponds to the trivial operating mode 
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of the heat flux splitter as a 50:50 splitter. 

On the other hand, if we change the Fermi level of 
disk 3 to Eps = 0.8 eV, we have an asymmetric situa¬ 
tion. The resulting powers V 12 and V 13 are plotted in 
Fig. 2 versus the radius R of the virtual sphere choos¬ 
ing a disk radius of D /2 = 100 nm. First of all, we see 
that V 12 can be almost three orders of magnitude larger 
than V 13 which means that nearly the full heat flux is 
directed from disk 1 towards disk 2. Here the heat flux 
splitter acts as a 99:1 beam splitter for distances around 
on micron. Clearly, by a gradual change of the Fermi 
level all values between 50:50 and 99:1 are attainable. 
Regarding the decay of the fluxes with R we find that 
in the far-held V oc 1/R 2 and in the near-held regime 
V oc 1/R 2 expected from the dipole-dipole interaction, 
where the held radiated by each disk decays as 1/R or 
1 /R 3 in far- or near-held regime, respectively. Of course, 
the direction of the heat how can be directed towards 
disk 3 by interchanging the Fermi levels of disk 2 and 3 
dynamically by electrical gating, for instance. 

To explain the huge difference of magnitude of the 
power echanged between two graphene nano-disks with 
identical Fermi levels and that between two disks having 
different Fermi levels, we analyze the spectra of the power 
huxes V 12 and V 13 calculated in Fig. 2. These spectra 
are plotted in Fig. 3(a) for the radius R = 2 D /with 
D = 200 nm, i.e. R = 231 nm. At small wavelengths, 
we see that both spectra are very similar. In contrast, 
at larger wavelengths, the shape of spectra are radically 
different. Excepted in a narrow spectral region cen¬ 
tered around the thermal wavelengths = ch/^ksTi) 
(i.e. spectral region where the Planck function is maxi¬ 
mal) , the spectrum of the exchanged power between two 
graphene nano-disks with different Fermi levels is more 
than two orders of magnitude smaller than that for two 
nano-disks with identical Fermi levels. Furthermore, the 
spectra show sharp resonances. These resonances are 
located exactly at wavelengths where the polarizability 
of the graphene disks become resonant (see Fig. 3(b)). 
These resonances are due to the localized plasmon polari- 
tons in graphene. At Ep = 0.1 eV this localized plasmon 
is located in the midinfrared at A = 25.1 fim. When the 
density of electron increases in graphene this resonance 
shifts to smaller wavelengths so that at Ep — 0.8 eV it 
is located in the infared at A = 7.75/im. When using 
identical Fermi levels these resonances overlap perfectly 
resulting in a large heat flux, while when using differ¬ 
ent Fermi levels this overlap becomes smaller resulting in 
smaller heat fluxes. Hence, we see that, by changing the 
Fermi level of graphene, we can switch the different chan¬ 
nels for the heat transport in this nanostructure network 
on and off. 

Finally, we show on Fig. 4 the angular dependence of 
heat flux exchanged between the disks for different sep¬ 
aration distances. We observe that the exchanged power 
V 12 decays monotically when 6 grows, whereas 7^13 in- 




Figure 3: (a) Spectral power exchanged between the graphene 
disks 1 and 2 (red) and between 1 and 3 (dashed blue) in 
near-held for the case considered in Fig. 2. (b) Modulus of 
the polarizability of the nano-disks from Eq. M for the two 
different Fermi levels Ef. 


creases monotonically. Moreover, when the mismatch of 
the Fermi levels in the graphene nano-disks is large the 
angular dependence is weak, whereas when the Fermi lev¬ 
els are identical the angular dependence becomes more 
important. 

In conclusion, we have introduced the basic build¬ 
ing blocks for a control of heat flow directions at the 
nanoscale by introducing a heat flux splitter based on a 
graphene nano-disk network. The working principle re¬ 
lies on the near-held interaction of the localized plasmons 
in graphene nano-disks which can be tuned dynamically 
by electrical means. In this Letter, we have limited our¬ 
self to the description of the simplest conhguration for 
a heat flux splitter. However, the basic concepts intro¬ 
duced here could be easily generalized to deal with heat 
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Figure 4: Angular depedence of heat power V 12 (red) and V 13 
(dashed blue) exchanged in the near-held between graphene 
disks of 100 nm radius in the heat flux splitter. The Fermi 
level of disks 1 , 2 and 3 are Ep 1 = 0.1 eV, Ef 2 = 0.1 eV and 
Ef 3 = 0.8 eV while their temperatures are T± — 350 K, T 2 = 
T 3 = 273 K. The radius of the virtual sphere is R = 500 nm. 

transfers in much more complex plasmonic networks. In 
such networks the near-held heat flux splitting could be 
used to dynamically manage heat exchanges in integrated 
architectures of nano-objects. 
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